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4^ Abstract 

Observations in Quantum Mechanics are subject to complex re- 
strictions arising from the principle of energy conservation. Deter- 
mining such restrictions, however, has been so far an elusive task, and 
only partial results are known. In this paper we discuss how con- 
straints on the energy spectrum of a measurement device translate 
into limitations on the measurements which we can effect on a target 
system with non-trivial energy operator. We provide efficient algo- 
rithms to characterize such limitations and we quantify them exactly 
when the target is a two-level quantum system. Our work thus identi- 
fies the boundaries between what is possible or impossible to measure, 
i.e., between what we can see or not, when energy conservation is at 
stake. 



1 Introduction 

The success of quantum computation [TJ [2] and quantum simulation [31 0] 
schemes depends in part on our ability to measure quantum systems with 
good enough precision. Sometimes (e.g., in ion-trap experiments [5]), such 
measurements are conducted in systems with a non-trivial energy opera- 
tor, and so are strongly limited by the law of energy conservation. Indeed, 
Wigner was among the first to notice the impossibility of measuring exactly 
any observable described by an operator which does not commute with the 
system's conserved quantities. He showed, nevertheless, that an arbitrarily 
close measurement of such observables was possible if the dimensions of the 
measurement device were large enough [6]. These results were formalized 
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in the Wigner-Araki-Yanase (WAY) theorem [7J |S] and quantified some time 
later by Ozawa [9], who provided a general uncertainty relation which bounds 
the mean square noise in the measurement of an arbitrary observable as a 
function of the variances of the conserved quantities in system and measure- 
ment apparatus. This relation has been applied successfully to estimate the 
error of certain lab-induced evolutions of two-level quantum systems [TUt [TTj. 
Sadly enough, and despite its generality and range of applicability, Ozawa's 
uncertainty relation is not tight, and sometimes greatly underestimates the 
errors it tries to bound. 

In this article we study how measurements of a target system effected by 
a given quantum device depend on the relation between the energy spectra 
of one system and the other under the assumption of energy conservation. In 
case the target is a two-level system, we find that the difference between the 
set of all conceivable two-level measurements and the set of all measurements 
implementable by measurement devices with a battery in one of the states 
B = {pb} can be quantified operationally by a simple parameter, t(B). We 
will calculate this parameter in two interesting scenarios: 1) measurement 
devices with finite energy spectrum, 2) measurement devices with unbounded 
energy spectrum, but average energy E. Additionally, we will provide an ef- 
ficient algorithm to characterize the set of accessible measurements for target 
systems of any dimension. Finally, we will propose a physical mechanism to 
explain why we should expect to estimate non-trivial observables when the 
spectra of target system and measurement apparatus are just approximately 
resonant. 

Before starting, though, we would like to call attention to the recent and 
related work of Ahmadi et al. [H], which also analyzes how the quantum 
state of the control device influences its capabilities for quantum measure- 
ments and control. Among its main results the reader can find a reformula- 
tion of the WAY theorem in the language of resource theories [T2l [T3] and 
the realization that, in finite dimensional systems, there does not exist such 
a thing as an optimal universal ancillary state, i.e., the optimal quantum 
state of the control device will depend on which specific quantum operation 
we wish to implement in the lab. 
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2 Measurements under energy conservation 



2.1 The measurement model 

A general measurement over a quantum system in state p G B(H) is described 
by a set of positive semidefinite operators {M x }™ =0 , with Ylx=o ^ x = ^ such 
that 

p(x) = ti(pM x ) (1) 

denotes the probability of obtaining the outcome x. In general, such a mea- 
surement is physically realized by attaching our target system S to a pointer 
P, and making them interact via a third system C that we will denote the 
clock. 

We will make the assumption that the clock that switches on and off 
the measurement interaction is a macroscopic object, and so the evolution it 
subjects systems S and P is approximately unitary. We will further assume 
that the clock does not transfer energy to the system under observation. Such 
is the case, for instance, when the interaction is introduced adiabatically. 
Indeed, let g(t)H]fp be the interaction effected by the clock on system SP at 
time t, with g(—oo) = g(oo) = 0. If g(t) varies slowly enough, the adiabatic 
theorem guarantees (modulo level-crossings) that the energy of system SP 
will be the same at the end of the process [15]. A subsequent von Neumann 
measurement of the state of the pointer will then return the outcome x of 
the effective measurement on system S. 

As we will see, energy conserving measurements are very limited, so we 
must include in the picture a fourth quantum system, the battery (5), to 
model the energy exchange between the target and the measurement appa- 
ratus. The battery has a non-trivial energy operator Hb, and only interacts 
with the system SP via the clock. 

The unitary interaction U between the target, the pointer and the bat- 
tery mediated by the clock will thus be such that the total energy H? of 
the system SPB will be conserved on average, i.e., {^\spbHt\^) spb = 
{iP\spbU^ HtU\iP) spb, for all joint states of target, pointer and battery. 
Now, the pointer can always be chosen such that its local energy operator is 
trivial- Denoting by Hs the energy operator of system S, this means that 

1 Indeed, given that the pointer must hold a reliable record of the outcome of the 
measurement, the states {|a)} must be orthogonal and do not evolve in time, i.e., they 
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H T = H S ® Ipb + hp ® H B . (2) 

For illustration, think of homodyne measurements in quantum optics: 
there, the target is the laser beam to be measured, and the battery is a high 
energy light pulse. The displacement of these two beams through optical 
fibers plays the role of the clock, which switches on and off an interaction 
with a beam splitter and two photodetectors. Finally, the difference between 
the intensity of both photodetectors is then printed on a piece of paper (the 
pointer). 

2.2 The role of the battery 

Suppose that we wish to measure our target system by means of a battery- 
less quantum device. Appendix |A] shows that the POVM elements {M x } 
describing any such measurement must commute with Hs- If, in addition, 
system S has a non-degenerate energy operator H s = J2 m E m \m){m\, then 
each M x must admit a diagonalization of the form M x = Yl m P^\ m )( m \- 
By the completeness relation, we further have that Y1 X P^ = 1 f° r an m - 
It follows that the measurement of any property x can be simulated by the 
following process: 

1. Measure the energy operator Hs, thus obtaining an outcome m (corre- 
sponding to the energy value E m ). 

2. Generate x randomly according to the distribution p(x) = ■ 

Measurements in this scenario are thus not very 'quantum', in the sense 
that all we need to return an outcome is the expression of the energy density 
of system S. In particular, no matter how much entangled system S is with 
another system S', two experimentalists at each site could never violate a 
Bell inequality or even prove that their joint state is entangled. Indeed, let 

must be eigenstates of the hamiltonian Hp of the pointer. If the energy of such states is 
equal, then we can assume that Hp = 0. If not, we can always attach a pointer P' to P 
with Hp/ = and make the joint system evolve via the unitary V = J2 a l a )( a l^ ® D(a) p> , 
where D(a)\k) = \k + a (mod n + 1)). If P' is initially set to the state |0), this unitary will 
copy the measurement information to P'. Moreover, [Hp <8> Ip<, V] = 0, i.e., V represents 
an energy-conserving operation. We can hence take P' to be our real pointer and regard 
P as part of the battery system. 



4 



Pss' be the joint state of systems S, S'. Then, if the experimentalist at S is 
restricted to perform measurements which commute with H$, it is easy to 
see that pss' will produce the same bipartite measurement statistics as the 
classical-quantum state 

Pss' = ^2p(m)\m}(m\ s <g> p% l) ■ (3) 

m 

where p(m) is the probability that system to has energy E m and p ( - m - ) = 
(m\p SS '\m)/p(m). 

How does the picture change when our measurement device has a battery? 
Let us assume that such is the case, and that H B = p, n \n)(n\ is non- 
degenerate. This last condition is not restrictive in the least, since we can 
always introduce degeneracy later via the pointer. 

As before, any effective measurement in the joint system SB will be 
described by a complete set of POVM elements {M x } commuting with the 
energy operator of the system. The difference stems in this case in that the 
new energy operator is not H s , but H s <g> l B + I5 <8> H B . 

This can make a huge difference. Coming back to the example of homo- 
dyne measurements, note that a measurement of the quadrature x cannot 
be performed via passive linear optical elements and photodetectors alone. 
Indeed, unlike in homodyne measurements, any such measurement would 
commute with the total photon number. Homodyne measurements can nev- 
ertheless be performed via linear optical elements up to arbitrarily good 
approximations if we introduce a second high energy laser pulse (i.e., a bat- 
tery). 

Note that 

H T = ^2{E m + p n )\m)(m\ <g> \n)(n\. (4) 

m,n 

We will say that there exists a resonance between the hamiltonians H B , Hs 
iff there exist to, to', n, n' such that 

E m — E m > = p n — p n i. (5) 

In the particular case that there are no resonances, the eigenspaces of 
Ht are given by {|to) <8> \n)}. Any operator M commuting with Ht will 
necessarily be of the form M = J2 m n M n>m \m)(m\ <8> \n)(n\. Let p B £ B(H B ) 
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be the initial state of the battery. Then, any measurement {M x } performed 
over the system SB through this scheme can be simulated by the process: 

1. Measure the energy of system S, obtaining outcome m. 

2. Output x randomly from the distribution p(x) = tr((m\M x \m) pb) ■ 

Non-resonant hamiltonians Hb thus do not provide any advantage towards 
measuring or interacting with system S in a quantum way. 

What about resonant hamiltonians? Suppose that our target system S 
is a two-level system, with Hs = |1)(1|, where |1) (|0)) denotes the excited 
(ground state) of S. Let us study which kind of effective measurements we 
can implement on system S if we couple it to a battery with hamiltonian 

HB = Ztlk\k){k\- 

The total energy of the system can be seen equal to 

h t = Etl Him ® \k)(k\ + 1 x >< 1 1 ® \k - i)(k -i\} + 

\d-l)(d-l\. (6) 
The eigenspaces of H T are thus defined by the projectors 

|0)(0|®|0)(0|, 

{|oxo|®|fcX*l + |iXi|®l*-iX*-i|}& 

\l)(l\®\d-l)(d-l\} r (7) 

Call the corresponding spaces % k , with %° = span{|0)|0)} and H d = 
span{|l)|rf — 1)}. Then, a generic POVM element in H k has the form 

M k = (M k ) ab \a}(b\®\k-a}(k-b\, (8) 

a,6=0,l 

where | — 1) = \d) — by definition. We will call such POVMs physical, since 
they correspond to the actual operation effected on the quantum system SB. 

Let pb be the state of the battery. Then, a physical POVM of the form 
M x = ® k M k induces in system S an effective POVM 

d 

{M x ) ab = J^ik - b\ PB \k - a}(M k ) ab . (9) 

fe=0 
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Now, take p B = \^b)(^b\, with \%j) B ) = Y^t=l \ k ) and let {Mjx be an 
arbitrary two-level POVM that we want to approximate. Set the physical 
POVM elements to be equal to {M x } x , i.e., (M^) ab = (M x ) ab . Then, one can 
verify that 

(M X U = {M x ) ab (l + ■ (10) 

It thus follows that any measurement can be approximated up to arbitrary 
precision by taking d large enough. It is straightforward to extend this result 
to multi-level target systems. 

Notice that, in order to perform non-trivial quantum measurements over 
the previous system, an exact resonance between Hb and Hs is required. 
This is certainly counter-intuitive: one would expect that ancillary systems 
with energy operator H B « Hb nearly resonant with Hs should induce 
similar effective measurements over system S (and thus approximate the set 
of all possible two-level measurements for d > 1). In Section [6] we provide 
a possible solution to this apparent paradox, by invoking the existence of 
hidden continuous degrees of freedom. 

The aim of the rest of the article is to determine exactly how the cardi- 
nality of the spectra and/or the energy of our measurement device constrain 
the set of effective POVMs that such a device is able to implement on its 
target system. But, before this, we will have to specify means to quantify 
such constraints. 



3 Quantifying the size of the set of accessible 
measurements 

Imagine that we hold a measurement device whose battery we can initially 
prepare in a set of states B = {pb}, and call Ai(B, d) the set of all effective 
POVMs which it allows to implement in a target system of dimension d and 
energy operator Hs G B(C d ). In the following section, we will try to quantify 
how Ai(B, d) compares with Ai(d), the set of all possible POVMs in C d . 

To do so, we must first introduce a natural distance between two differ- 
ent POVMs. Suppose that we have a device capable of implementing either 
the measurement M° = {M^} x or M 1 = {M]} x . with probability 1/2. We 
let it act over a suitably prepared quantum state p and, from the outcome 
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x obtained, we try to guess which of the two POVMs our machine is actu- 
ally implementing. It can be shown (see Appendix [U]) that the maximum 
probability Pq of correctly guessing the POVM is then given by 

P c = i{l + dist c (M°,M 1 )}, (11) 
where distc(M°, M 1 ) is defined as 

dist c (M°, M 1 ) = ~ max]T |tr{p(M x ° - M*)}|, (12) 

X 

and the maximization is to be performed over all normalized quantum states 
p. We will call this latter expression the classical distance between POVMs 
M°, M 1 . It satisfies the triangle inequality (i.e., it is a proper distance), and 
has maximum value 1. The reason why we call it classical is that, in the 
previous protocol, the POVM is guessed by analyzing the classical data x. 
The classical distance is somehow related to the distance between quantum 
maps introduced in |16j . 

Analogously, we can define a quantum distance between POVMs, by 
means of a protocol where the player is allowed to input part of an entangled 
state pDQ in the measurement device, which then performs a demolition mea- 
surement on system D. Depending on the outcome x of such a measurement, 
the player will implement a POVM N%, with outcomes a G {0, 1} on system 
Q in order to decide which of the two POVMs, M° or M 1 , is actually being 
measured. The probability Pq of correctly guessing the POVM can then be 
seen equal to 

P Q = ^{l + dist (3 (M ,M 1 )}, (13) 

with 

dist Q (M°,M 1 ) = imaxV||tr D (p DQ (M°-M ;E 1 )®lQ|| 1 (14) 

I PDQ Z ' 

X 

where pdq varies over all possible states pdq € B(T-Lr> ® Hq) and all pos- 
sible Hilbert spaces %q. See Appendix [C] for a proof. As the classical 
distance, distg satisfies the triangle inequality and has maximum value 1. 
Also, dist Q (M°,M 1 ) > dist c (M°, M 1 ), for all M°, M 1 . dist Q (M°,M 1 ) ac- 
tually corresponds to the diamond norm [T7] between the quantum channels 
fia (p) = £x tr G° M *Mz|> wi t h a = 1,2. 
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It can be seen that, if M°, M 1 are two-outcome POVMs, dist Q (M°, M 1 ) = 
dist c (M°,M 1 ). However, even in dimension 2 there exist examples of POVMs 
where dist Q (M°,M 1 ) > dist^M^M 1 ) (Appendix EJ). 

These two distances suggest a simple way to quantify the difference be- 
tween a particular set of POVMs M.' acting on C d and the set Ai(d) of 
all possible quantum measurements in that space, by computing the maxi- 
mum distance between an arbitrary element of Ai(d) and the set M.' . This 
intuition leads to the following definitions: 



e c (M')= max dist c (M, M'); en(M') = max distn(M, M'). (15) 

MeM(d) M&M(d) 

ec(M.'), eg(A'f) thus measure the worst-case probability of correctly dis- 
tinguishing a general POVM M e Ai(d) from its closest element in Ai' in 
classical and quantum protocols, respectively. 

We have just defined two quantities to measure the performance of a quan- 
tum measurement device. Our next step will be to evaluate such quantifiers 
when the target system is a two-level quantum system. 



4 Two-level systems 
4.1 General considerations 

As in section [2J suppose that our target system S is a two-level system with 
Hs = A|l)(l|, which we couple to a battery B with hamiltonian Hb- 

For the time being, assume that the spectrum of Hb is discrete. We say 
that an increasing sequence of eigenvalues of Hb forms a chain of length L iff 
it can be written as (z/ + A/c)^~q, for some v G 1R. We call a chain maximal if 
it is not a subset of a larger chain. It is clear that each /z 6 spec (_£/#) belongs 
to a unique maximal chain. Hence we have that 



Ht = £^|0)<0|®|j,0)(j,0| + 

+ Efc2~ V; + fcA){|0X0| <8> \j, k)(j, k\ + |1)(1| 9 \j, k - k-l\} + 

+{uj + L(j)A)|l)(l| ® \j,L(j) - l)(j,L(j) - 1|, (16) 

where each j denotes a maximal chain of length L(j) and \j, k) is the nor- 
malized eigenvector of Hb with eigenvalue Vj + kA. 
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The eigenspaces of H T are thus 



|0>(0|® |j,0)(j,0|, 
{|0)(0| ® \j, k)(j, k\ + |1)(1| ® |j, fc - l)(j, k - 

(17) 

In analogy with the previous section, call these spaces W' , with "H J '° = 
span(|0)|j, 0)) and W' L ^ = span(|l)|j, L(j) — 1)), and let ps be the state of 
the battery. Then, a physical POVM {M x = ® k ,j M x' j }, with 

M** = J2 ( M ^U\a)(b\®\j,k-a}(j,k-b\, (18) 

o,6=0,l 

induces in system S an effective POVM 

(M x ) ab = J2H WpbIj, k - a) (j, k - b\}(M^) ab . (19) 

j k=0 

If H B is of the form (T2"U|) . by eq. (ITU]), a measurement device possessing a 
battery B' with energy operator Hb 1 = Yln=o n \ n )( n \ (where L = maxj L(j)) 
can simulate the above measurement via the following protocol: 

1. Choose j randomly according to the distribution {pj}, with pj = 

2. Prepare the state pj oc Ylfki^odi ^IpbU, k')\k){k'\ in the ancillary sys- 
tem B'. 

3. Measure the joint system SB' with the physical POVM {Ml = ® k M^ k } 
(after any arbitrary completion, note that ^Mjj. < I). 

4.2 How close is our measuring device from being per- 
fect? 

Imagine that our preparation devices are capable of setting up any battery 
state pb G B at the beginning of the experiment. Suppose further that the 
energy operator of system B is of the form 
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i(3)-l 

Hb = J2 E fa + *A)|*,jXM- (20) 

j fc=0 

In Appendix [D] we show that, given those conditions, 



e c (M B ) = e Q (M B ) = ^{l-r(B)}, (21) 



where 



r(S)=max^ £ fc + l\p\j, k)\ . (22) 

p j fc=0 

Moreover, let p* G £> be any state maximizing ( 1221) . Then, for any general 
one-qubit POVM M = {M x } x £ M, p* allows to generate an effective two- 
level POVM M = {MJ, with 

(M x ) ab = (M x ) ab {[l - r(B)}5 ab + t(B)}. (23) 

The resolution of ( |22i) thus allows to define a prepare-and measure strategy 
which simulates the whole set Ai with accuracy 1 — t{&). 

Correspondingly, if the ancillary system has a continuous energy spec- 
trum, the value of r(p B ) equals: 

t( Pb ) = J dE(E\p B \E + A). (24) 

In the particular case that p B is a pure state with energy density p(E)dE, 
we can re-express this last equation as 

r(p B ) = J dEp x l\E)p x l\E + A). (25) 



-E 



For example, take p(E)dE w s ^=dE. Then, r(ps,A) ~ e ; i.e., the 
measurement device state is useless to measure systems with A 3> *Jo ' . 
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4.3 Finite spectrum 

Picture an experimental scenario where we only have control over the first 
d energy levels of our battery system. That is, we can only prepare battery 
states with energy density p{E) = Ylk=oPk$(.E ~ Ek), for some energy values 
{E k }. Equivalently, suppose that the spectrum of H B is finite, i.e., that 
any measurement of the energy of the battery can only produce a finite set 
of outcomes. Denote by d the cardinality of spec(ife). What is the set 
of measurements Ai(C d , 2) which one can perform over system S with this 
model, and how does it differ from the set of all possible POVMs M.(2)7 Most 
importantly, how fast does Ai(C d , 2) tend to Ai(2) in the limit d — > oo? 

As we saw in section 14.11 the existence of more than one chain in the 
spectrum of the battery does not provide any advantage. In the following we 
will hence assume that H B = Y,£*kA\k)(k\. 

Our aim is therefore to study the set of effective POVMs implementable 
via measurement apparatuses equipped with a battery of energy distribution 
{Pk > : Ylk=oPk = !}• By Appendix IB"1 we can take the battery states 
to be pure, i.e., = Ylt=o Ck \^)- Then, according to the last section, in 
order to quantify the difference between JH(C d ) and Ai, all we have to do is 
compute 

{d-2 d-1 "| 

^|c fc c fc+1 | :^| Cfc | 2 = l I. (26) 
k=0 k=0 ) 

In Appendix lEl we show that the solution of (1261) is cos (5x5;) • It follows 
that 

e c (M(C d )) = e Q (M(C d )) = \{l- cos (j^) } . (27) 

Curiously enough, the states \ip* d ) maximizing eq. ff22|) happen to have a 
non-trivial energy distribution, given by 

p(k) = - A- sin 2 ( ,k = 0, d-1. (28) 



d+1 \ d 
4.4 Finite energy 

Consider now a scenario where in principle we can prepare any initial battery 
state, but we do not wish to spend too much energy in the process. Note 
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that the energy of states |^) in the previous section is equal to (d — l)A/2, 
i.e., it grows linearly with the dimension d. This makes one wonder whether 
such energetic states are actually necessary in order to attain a good approx- 
imation to Ai(2). Or, in other words, how well can we approximate Ai(2) 
when our battery is infinite dimensional, but its energy is bounded? 

Call M(E, 2) the set of two-level POVMs attainable via devices with a 
battery of energy smaller or equal than E > for Hb = Ylk=o kA\k)(k\. In 
analogy with the previous section, define e^(E), €q(E) as ec({p : tr(p_£/s) < 
E}), e<g({p : tr (pHs) < E}), respectively. From the previous section and 
the observation that (^l-ffsl^) = (d— l)/2, it follows that M(E, 2)e_> 00 — > 
M{2). Moreover, 

e* c (E) = e* Q (E) <^1- cos J — j } . (29) 



that is, for large E, Gq(E) < 0-6168A 2 ^y e wonc [ er w hether this last bound 
is tight, or close to tight. Could it be that there exists an energy threshold 
E > such that M(E, 2) = M(2) for E > E ? 

Again, by eq. (122]) . the answer to these questions depends on how well 
we can approximate 1 by X^fcLo l c fc c £:+i|; under the constraints J2 k \ck\ 2 = 1, 
J2 k \c k \ 2 kA<E. 

What we find in Appendix [F] is that r({p : tv(pH B ) < E}) = <p(E/A), 
with 

/ \ . z + {fi: = 2A} 
ifi(z) = mm ^ , (30) 

where j n> i denotes the first zero of the Bessel function of the first kind J n . 
The numerical minimization of this function is not problematic, since all its 
extreme points are global minima, see Appendix |F] For z 3> 1, f(z) behaves 
as 

0.9468 

~ 1 ( 31 ) 

Consequently, 

1 - 4734A 2 

eh(E) = e* c (E) = -{1 - V (E/A)} « — for E » A. (32) 
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The states |^^) = YlT=o c k\k) with average energy E maximizing eq. 

will be denoted power states. Their coefficients {c^} are given by the 
recurrence formula: 

k + {/i : = 2A*} 

Cfc+l — ^ Cfc — Cfc_i, (00) 

where A* G 1R + denotes the minimizer of (1501) . 

It is instructive to compare how useful power states are for high precision 
measurements as opposed to the usual coherent states, of the form \a) = 

e~4- Y^h=o vfyl^)' with E = |a| 2 A. For coherent states, 



00 Ul 2fc+1 1 A 

r(| a )) = V ' 1 ^l-— 1 — = 1- it, (34) 

where the approximation holds for |a| 3> 1 [IE]- Now, we can quantify 
measurement precision by log 10 (l — r), i.e., by the number of significant 
decimal places up to which we can approximate M.{2). According to this 
definition, the last equation means that, for E/A ^> 1, power states double 
the precision we may reach with coherent states of the same energy. It is 
hence an interesting question whether current technology allows producing 
power states in the lab. 



5 Characterization of the sets of accessible 
measurements 

The results of the last section show that M(C d , 2), M(E, 2) ^ M(2), and 
determine an operational distance between these sets. They leave open, 
though, the problem of characterizing which two-level POVMs M 6 M(2) 
can be realized with a measurement device with bounded energy spectrum, 
or average energy. More generally, given a POVM M e M(d), we would 
like to decide if such a POVM can be implemented in a d- dimensional target 
system S with non-trivial energy operator Hs with the aid of measurement 
devices with energy operator Hb- 

The dual of this problem is also of interest for the quantum information 
theory community: given an arbitrary vector of operators V = (Vo, V\, V n ), 
maximize 5^" =0 tr(M x V^), for M 6 A^C^). This class of problems includes 
the energy constrained state discrimination problem: given a device which 
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randomly prepares a multi-level quantum state {pj}™ =0 according to some 
known probability distribution {jOj}™ =0 , maximize the probability of correctly 
guessing which state pi was produced by measuring it with a device possessing 
a ci-level battery. 

In Appendix[B]we show how to formulate the characterization of Ai (C d , 2) 
and linear optimizations over Ai(C d , 2) as semidefmite programs [TH] involv- 
ing 0(d) 2x2 complex matrices. Due to this small scaling, using standard 
convex optimization packages like sedumi [21], we found that a normal desk- 
top can carry optimizations over A4(C d ) for d > 4000. We used such pro- 
grams to prove that the state p* with energy distribution ( )28|) is not a univer- 
sal resource state, in the sense that there exist measurements M G A4(C , 2) 
which are not achievable with such an ancillary state, thus recovering the 
conclusions of [H]. As shown in Appendix [Bj the algorithm can be easily 
adapted to characterize the set A4(C d , d') for arbitrary (given) hamiltonians 
Hb G B(C d ), Hs G B(<C d ). The corresponding semi definite program involves 
0(dd') d' x d' complex matrices, and is therefore computationally efficient in 
both the ancilla and the target dimensions. 

The next problem is how to characterize Ai(E, 2). Since a full charac- 
terization of Ai(E,2) would involve optimizations over infinite dimensional 
energy distributions {pk : k > 0}, it is unlikely that we can reduce it to 
a semidefmite program. The approach followed in Appendix [B] is to define 
two sequences of inner ({M d (E, 2)}) and outer ({M d (E,2)}) approxima- 
tions of M(E,2), i.e., M d (E,2) c M{E,2) c M%(E,2) c M{2), with 
lim^oo M^E^) = lim^oo M%(E,2) = M(E,2). Each of these approxi- 
mations can be computed via a semidefmite program involving 0(d) 2x2 
positive semidefmite matrices. Moreover, the speed of convergence of the 
scheme is bounded by 0(E/dA) (see Appendix [B]). Maximizing a linear ex- 
pression f(M) over M e Ai(E,2) can thus be accomplished by optimizing 
over the sets A4®(E, 2) and Ai ! d (E, 2). The first optimization will return an 
upper bound f° on the maximal value / = max M& M(E,2) f(M), while the 
second optimization will return a feasible POVM M G Ai(E, 2) attaining 
a sub-optimal value fj. The difference f d — f\ will quantify the error of 
the approximation. This algorithm can be adapted as well to describe ef- 
fective measurements in multi-level quantum systems with devices of energy 
bounded by E if the spectrum of Hb is either finite dimensional or admits a 
sufficiently simple description. 
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6 Why can we make quantum measurements 
at all? 



As we pointed out in section [21 in order to perform non-trivial quantum mea- 
surements over system S, the hamiltonian of our device's battery Hb must 
be exactly resonant with the hamiltonian of the system under observation 
Hs = A|l)(l|. Such a state of affairs seems very unphysical, since it would 
allow us to distinguish a hamiltonian A|l)(l| from, say, Hb = (A + e)|l)(l| 
for arbitrarily small e ^ 0. How can we solve this paradox? 

One possibility is that the energy difference is accounted for by extra 
degrees of freedom with continuous energy operator. Suppose that we are 
actually trying to perform a non-trivial measurement over a two-level sys- 
tem S with Hs = A|l)(l| by coupling it to a two- level battery with energy 
operator H B = (A + e)|l)(l|, with A > e > 0. If \ip B ) = c |0) + ci[l) and 
e = 0, we would expect to find a t(iPb) — l c o c i|- F° r e > 0, on the con- 
trary, the theory predicts r(ip B ) = 0. Imagine, though, that there is an extra 
degree of freedom in the lab, i.e., a third quantum system C with wavefunc- 
tion \ifj c ) = jo°° v 1/2 (E)dE\E), where a(E)dE denotes the energy density of 
system C, and {\E)dE} is an (undercomplete) orthonormal basis of energy 
states. This continuous degree of freedom could correspond, for example, to 
the kinetic energy of the clock mediating the measurement interaction. As 
shown in section H~2| the value of r for a pure ancillary state \4>b) with energy 
density p(E)dE used to control a two-level system with energy operator Hs 
is given by eq. fl25|) . 

The energy density p(E)dE of the system BC in the previous case can 
be seen equal to 

p(E) = {cr(E)\c \ 2 + a(E — A — e)\ Cl \ 2 }dE, (35) 

and, consequently, the potential for measuring system S is determined by 
the quantity 

t{\^b) ® m) = J dE(a(E)\c \ 2 + a(E - A - e)| Cl | 2 ) 1/2 • 

■{a{E + A)|c | 2 + a(E - e)\ Cl \ 2 ) l l 2 dE. (36) 

In the particular case in which a(E)dE is very concentrated around the 
energy values [E — a, E + a], with e a A, the above expression is 
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very close to | CqC^ I , as expected. Note, though, that due to the localization 
condition, r({ipc}) ~ 0, i.e., the auxiliary system C cannot do much by itself. 

As an example, imagine that \iPb),Ps correspond to the quantum descrip- 
tion of the inner degrees of freedom of two systems, say, molecules, which 
are able to move in space. Suppose further that such systems are, at the be- 
ginning of the experiment, in a very- well localized state in momentum space 
(perhaps because they have been sent to collide and thus interact). Call 
|V , )c, \^)d the corresponding wavefunctions, and assume, for simplicity that 
they are identical and equal to 

\<j>) oc J d 3 pe~£\p). (37) 

Under the assumption that the masses of each system are both equal to m, 
we have that 

[ if p 2 - 

\i>c) ® |<M oc / dP T dpe-^e-^\P T )\p), (38) 

where pc = \Pt + P, Pd = \Pt ~ P- If we further assume conservation of 
total linear momentum, we can imagine a prior measurement of Pt followed 
by an interaction between the systems A, S and the remaining of CD. Now, 
for whatever outcome of Pt, the energy density of the remaining of system 
CD (we only count the kinetic terms) is given by the distribution 

a{E)dE=( — ) —e~ — E 1/2 dE. (39) 
V cr / y/n 

The application of formula (l36l) gives us 



t(ABCD) >e~^ /^f^) 3/2 A| Co || Cl | e -^ E i/4 (jE ;_ e) i/4 

> | c 1 1 ci | e 2- . (40) 

The first inequality is close to being saturated for ^ C A, the second one is 
just an approximation. In either case, one can see that, as e — > 0, we recover 
the exact result in a continuous way. 
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7 Conclusion 



In this paper, we have investigated how the energy spectrum of a quantum 
measurement device limits what can be observed in a target system S. More 
concretely, we have quantified exactly the efficiency of measurement devices 
acting over a two-level target system as a function of their energy content and 
dimensionality. We have also provided efficient algorithms to characterize the 
set of effective measurements given the hamiltonians Hs, Hb of the target and 
the battery. Finally, we have identified a mechanism that allows the device 
to perform non-trivial measurements when the hamiltonians H B , H s are not 
exactly resonant, by invoking the influence of continuous hidden degrees of 
freedom. 

In future communications we will address the problem of quantifying the 
efficiency of measurement devices on multilevel systems, and characterizing 
the set of possible evolutions on system S induced by quantum devices with 
battery B. 
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A General remarks about energy conserva- 
tion 



A.l State evolution under energy conservation 

Let 1 and 2 be two quantum systems, described through Hilbert spaces 
%i,%2; respectively. Call H G B(%) the energy operator of system 1, and 
assume that the energy operator of system 2 is trivial (i.e., a multiple of the 
identity). We will next derive the set of possible transformations which we 
can induce in system 1 by coupling it to system 2 under the assumption that 
energy is conserved. Suppose then that we let both systems evolve according 
to a joint unitary operator Ui 2 , i.e., the system undergoes a transformation 
of the form p 12 — > Ui 2 pi 2 U\ 2 . 

The assumption that energy is conserved translates as 

tx{U 12 p 12 U\ 2 H ® I 2 ) = tr {p l2 H g> I 2 ), (41) 

where pi 2 G B(Hi <S> Hi) is the initial state of the bipartite system. Since 
this condition applies to any initial state p i2 , we have that U\ 2 (H <g> l 2 )Ui 2 = 
H <S> I2, or equivalently, 

[U, H®I 2 ] = 0. (42) 

If the initial system is in a product state, i.e., p\i = p\® IV^XV^I (w.l.o.g. 
we take the auxiliary system to be in a pure state), then, upon tracing out 
system 2, system 1 will have undergone an evolution pi — > ^2 i A i piA\, with 
{Ai} = (Ii g) (i^f/^Ii <g) l^)- Here denotes any orthonormal basis of 

7{ 2 • Note that 

[Ai, H] = [A\, H] = 0, Mi = h- (43) 

i 

Conversely, for any set of operators {Ai} C B(Hi) satisfying (T4"3"j) . there 
exists an ancillary system 2, a state |-?/> 2 ) G H 2 and a unitary Z7 12 G ® 
7£ 2 ) with [C/ 12 ,F«)l2] = 0, Ei^iPi4 = tr 2 (C/i 2 pi® IM^I^)- Indeed, let 
H = ^2 n E n P n , with {P n } being a complete set of orthogonal projectors. If 
[Ai, H] = [A\,H] = 0, then A; = © n A-"\ with supp(A (n) ) C supp(P n ). Call 
1-6^ the support of P n , and consider the operator U\ 2 = (B n U^ n \ with 
being unitary operators defined in the subspace %^ ® H 2 as 
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f/ 1 ( ?l^)i|0) 2 = ^4 (n) |^)iK)2, (44) 

i 

for all \ip n ) e H {n) . Obviously, [U 12 ,H] = 0. Moreover, it can be verified that 
Z i A iPl A\ = tT 2 (U 12Pl ®\0)(0\Ul 2 ). 

A. 2 Measurements under energy conservation 

As before, consider two quantum systems: system 1, our target (with energy 
operator H e B{T-L\)) and system 2, the pointer (with trivial energy opera- 
tor). We will model the measurement process of the system 1 via coupling it 
with system 2 through a unitary U\ 2 , followed by a von Neumann measure- 
ment of the pointer. As we saw before, after a fine-grained von Neumann 

measurement { |*)(* | }f=o 1 °f system 2, the state of system 1 will be ^— L , 

with probability pi = tr(pA\Ai). Here Ai = (Ii <S> {i^UuIi <S> \ip2), and so 
[A,H] = [Al,H] = 0. 

The statistics of a general von Neumann measurement of system 2 with 
outcomes x — 0, ...,n can be simulated by identifying disjoint subsets Q(x) 
of outcomes % = 0, d 2 — 1 with the effective outcome x. The probability of 
obtaining outcome x in any such measurement is given by 

tT(pA\A i )=ti(pM x ), (45) 

ieci(x) 

where the operators M x = J2ien(x) ^\Ai satisfy 

M X >0,J2M X = J2 A l A i = J i> t M - # 1 = °- (46) 

z i 

Conversely, any measurement associated to a set of operators {M x }" =0 
satisfying the above conditions can be realized via an energy conserving uni- 
tary interaction with a pointer followed by a von Neumann measurement 
of the latter. To see why, define A x = \[M X and apply the results of the 
previous section. 
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B Algorithms 



In this Appendix we will devise a collection of algorithms to characterize 
the effective POVMs generated by measurement devices with different con- 
straints in their batteries's energy distribution. 

Note that, for a fixed battery state ps (not just pure) eq. (1T91) can be used 
straightforwardly to formulate a semi definite program [19J to characterize the 
set of effective POVMs generatable via ps- However, if we allow ps to vary 
over a continuum of possible states B C B{%b), eq. f|T9|) is no longer an 
SDP. If we want to optimize over non-trivial sets of battery states, it is thus 
necessary to find a more appropriate reformulation of the problem. 

Given an energy distribution {pj )k > : ^ - Yuk=o Ps',* = 1}> consider 
first the problem of characterizing all one-qubit POVM elements which can 
generated by any possible battery with such an energy distribution. 

First of all, we can assume the battery state to be pure. Indeed, suppose 
that pb is mixed, and consider a purification \i/j)bb' in an extended space 
Hb ®Wb'- Then, we can write \^)b as 

L(j)-1 

\^)b = Y1 VHk\fJ,k), (47) 

3 k=0 

with \(pj t k) being eigenvectors of the energy operator Hb <8> Ib'- By redefin- 
ing the physical POVMs as M£ k (\(pj, k )(<Pj,k\ + \<Pj,k-i)(<Pj,k-i\)M% k g) 
^■B>{\ ( Pj,k)( ( Pj,k\ + l^j.fe-iX^'.fc-il), w e end up with a pure state \iP)b, with the 
same energy distribution as ps, defined over a battery system with a non- 
degenerate hamiltonian H' B with the same spectrum as Hb, and which allows 
to reproduce the same efffective POVM {M x } x . 

Now, the effective measurements generated by apparatuses with battery 
state \i[))b = J2j,k c j,k\j, k) are given by 

M, = EE c lk- a Cj, k- b {Ml> k ) ah , (48) 

j k=0 

where Cj-i = CjUj) = 0. 

Note that we can rewrite this equation as 

3 k 
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with 



k - k s ( t c-L ) «** ( ? ,„!., ) • < m > 

What conditions must {M^} satisfy? By definition, they are positive 
semidefinite. Also, by the completeness relation, they must satisfy 

where {pj t k = l c j,fc| 2 } denotes the energy distribution of the battery. 

Conversely, for any set of non-negative operators {M ] x ,k } x satisfying fl51~j) . 
one can find a physical POVM {M^ k } x such that Mp satisfies eq. fl50|) . 
A characterization of the set of POVMs attainable with states of energy 
distribution {pj,k} (remember that = Pj,L ( j) — 0) is thus given by 

m 

j k=0 

s.t. Mi* > 0,^ = ( P ^ ^ ) , (52) 

which, in finite dimensions, constitutes a semidefinite program (SDP) |19j . 

In finite dimensions, it is thus immediate to design an algorithm to char- 
acterize all effective POVMs generatable through ancillary states with en- 
ergy distribution defined via the shadow of a spectrahedron, i.e., all col- 
lections of numbers {pj t k} such that fc A?,feP.j,fc + Yli^iVi — f° r some 
matrices Aj^, Bi and some extra (free and fixed) variables y. Simply turn 
{pj,k} into free variables in program eq. (152]) and add the extra constraint 

£,-,* A j,kPj,k + B m > o. 

B.l Characterizarion of M(C d , 2) 

As remarked in section 14. 3\ we just have to consider batteries with energy 
operator Hb = Ylt=ok\k){k\. The energy distribution of the states of the 
battery {pk}t=o is j us ^ limited by the constraint £ fc IoPfc = 1. A program to 
characterize Ai(C d , 2) is thus 
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d 

M x = Yl u ' • 

fc=0 

s.t. M*>0,^M* = 

X 

d-l 

P-l =Pd = 0,^Pk 
fc=0 

B.2 Characterization of M{E, 2) 
Our aim is to solve the feasibility problem 

oo 

M x = Yl V/ '- 

fc=0 

,, «*>o, ), 

oo oo 

p-i = 0, 5^p fc = 1, ^ p fc A;A < £. (54) 

fc=0 fc=0 

Unfortunately, the above problem is not an SDP, since it involves an infinite 
number of free variables. 

The way we will solve this problem will involve relaxing or restricting the 
above program via truncations of order d in order to get inner and outer 
approximations to the set Ai(E,2). We will then show that such approxi- 
mations converge to Ai(E, 2) and derive some bounds on the speed of con- 
vergence. 

First, consider the obvious restriction M d (E, 2) = M(E, 2) n M(C d , 2). 
This is equivalent to program ([55]) with the additional restriction Ylk=oP k ^^ — 
E. This program thus constitutes an inner approximation of Ai(E, 2). 

As for a relaxation of (1541) . fix a number d G N greater than 1, and note 
that, for any energy distribution {pk} we have that 

kd * ^0 Vk-i J V P d + Pd ^ J 
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Pk 
Pfc-i 



= 1. 



(53) 



with P d = l- J2k=oPk- Also > notice that Y^k=o k Pk + dp <i < Y^kLo k Pk < f • 
A suitable relaxation of (1541) is thus the program 

d-l 

M x = ^Tm^ + M^, 

fc=0 

s.t. M k x > 0,fc = 0,...,d 

^ ^ Pd + Pd-i / 

P-i = 0,^p fc + P d = l,^)A;p fc + dP d <^. (56) 

fc=0 fc=0 

Let us call M°(E,2) D M(E, 2) the outer approximation of A4(i?, 2) 
so generated. An immediate question is how close M.®(E, 2) and Ai^E, 2) 
actually are. It turns out that, for any POVM M 6 M.°{E, 2) it is fairly 
easy to construct an approximate POVM iV 6 A4 d+1 (-E, 2). Indeed, define 
the matrices: 

M d = p ^ ),M, d = 0forx^0, 

^° +1= (o pj>^ +1 =0forx^0. (57) 
Then the POVM given by 

d-l 

N x = Y,M k x +M d x +M d x +1 (58) 

fc=0 

is trivially in M{E, 2) n M(C d+1 , 2). 

Given an arbitrary POVM M e M^(E, 2), its distance with respect to 
M T d+ i(E, 2) is bounded by 

dist Q (M,A^ +1 (£)) =i max _ V ||tr D {p D0 (M a! - JV S ) <g> I Q }||i 

2 PDQ ,NeM' d+1 (E)^ 



[ P * ° ),V*<d-l, 
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< - max 

2 pbq 



Mi - M; 



< 



"£\\tr D {p DQ (M d+1 

x 

\ max tr { PDQ (M d x +1 + M d x + M d+1 



< - 



2 PBQ 
1 



2P d 








< 



E 



A(d-l) 



(59) 



where the last inequality comes from the fact that (pd-i + Pd){d— 1) < E/A. 

Being this a rough estimation of the distance between Ai®(E,2) and 
A4^ +1 (E,2), we recommend the reader to run both programs in order to 
perform linear optimizations over Ai(E, 2). 



B.3 Extension to multi- level systems 

The above results can be extended to characterize the set of effective POVMs 
implementable in a target system of dimension d! and (non-degenerate) en- 
ergy operator Hs by a measurement device with (non-degenerate) hamilto- 
nian Hb- Indeed, given Hs = E m \n)(n\ and Hb = ^2 c Hn\n)(n\, the total 
energy operator is described by eq. (j3J). Re-expressed in terms of energy 
subspaces, the operator reads: 

H T = Y,E S P S , (60) 

s 

where 

d(s)-l 

P s = 2j \m(s,i))(m(s,i)\®\n(s,i)){n(s,i)\. (61) 

Here, d(s) < d! denotes the rank of the energy subspace P s , and {m(s, i)}f=Q 1 
{n(s, i)}i=o 1 are used to denote the eigenvectors \m(i,s)) <8> \n(i,s)) of Ht 
with eigenvalue E s . Note that (m(s, i)\m(s, j)) = (n(s,i)\n(s,j)) = 0, fori 7^ 
j. Now, let p B = |*bX*b|, with |* B ) = J2 S EtLT' c n(s,i)\m{i, s))® \n{i, s)). 
Any physical POVM implemented over the joint system SB must be of the 
form {M x = Q) S M X } X , with ^2 X M X = P s . Then one can check that the 
effective POVM performed on system S is given by 

M. = X)A& (62) 
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with 



( M £)m(»AM'j) = C *s,i C n(s,j)(M°) m (s,i),n(s,j)- ( 63 ) 

As before, it follows that {M*} x are simply limited by the positive semidef- 
initeness condition and the completeness relation 



= diag(p n(S)i) ), (64) 



X 



where p n (s,i) = \c n (s,i)\ 2 denote the energy occupation numbers of p&. 

Characterizing the set of all measurements effected on 5" by measurement 
devices with energy occupation numbers describable as the shadow of a spec- 
trahedron can then be trivially formulated as a semidefmite program. That 
includes the case where such numbers are just constrained by summing up 
to one, i.e., we can easily characterize the set Ai(C d , d'). The case where the 
spectrum of the battery system has infinite cardinality can be also attacked 
by hierarchies of SDPs as in the previous section. 

C Interpretation of dist^, distg 

Imagine a quantum device capable of performing a certain demolition mea- 
surement M a e Ai, thus producing an output x. We are further promised 
that the device is actually measuring either M° e M or M 1 e Ai with prob- 
ability 1/2. Suppose, then, that we input a given quantum state p and obtain 
an outcome x. Let p a p {x) = tr{pM°}. Clearly, the strategy which maximizes 
the probability of guessing a is to choose a = argmax{^(x) : a = 0, 1}. The 
maximum probability of guessing a by classical means is thus 




X 




X 




(65) 



with distc defined as in eq. (jl2]h 
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Suppose now that, rather than analyzing the inputs of the device, we let 
it measure subsystem D of an entangled state pdq E B{1-L d ®'Hq), where %d 
(jHq) denotes the Hilbert space corresponding to the device (the auxiliary 
system). After obtaining output x from the machine, we perform a POVM 
N x = {N® : a = 0, 1} over system Q, whose outcome will be our guess on 
the value of a. The maximum probability of success of this scheme is 



P Q =o m 0f E tr M M ° ® N ° + Ml ® N? )} = 

Z p,N x ' 
x 

L A P,N X * — ' 
x 

= \ + \ max E tr +^° -pI)i ( 66 ) 

with p a x = tro(pM£ ®1q) and tr + (T) denoting the sum of the positive eigen- 
values of operator T. Now, it is easy to see that, for any self-adjoint operator 
T, ||T||i = 2tr + (T) — tr(T). Substituting in (|66|) . we have that 



X X 

= i{l + dist Q (M°,M 1 )}, (67) 

where distq is defined as in CHI) and in the last line we made use of the fact 

that Zx ^Q(PI) = tr (P Ex M x ® h) = h for a = 0, 1. 

A natural question is whether distc and distQ are actually different. They 
can be proven equal if M°, M 1 happen to have only two outcomes. Indeed, 
suppose that such is the case. Then we have that 

dist Q (M°,M 1 ) = max \\tr D {p DQ (M° - M^jh- (68) 

PDQ 

Using standard identities of the trace norm [22], we have that 



\\tT D {p DQ (M° -Mo 1 )}!! = max tr{p DQ (M ° - M 1 ) ® A} 



KA<I 



< HAfJ-AfJlU. (69) 
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On the other hand, we have that 

dist c (M°, M 1 ) = max |tr{p(Af ° - M%)}\ = ||M ° - MH^. (70) 
p 

Comparison of these two relations yields distg(M°,M 1 ) = distc(M°, M 1 ). 

However, even in dimension 2 we can find pairs of POVMs M°, M 1 with 
dist Q (M°,M 1 ) > dist c (M°,M 1 ). Take, for instance the continuous POVMs 
defined by 

Mj = 2|^|#,M; = I 2 ^, (71) 

where dip denotes the invariant measure on the pure states of C 2 . Since both 
POVMs are invariant under rotations, we can choose p = |0)(0| in eq. ffl2|) . 
Then we have 



dist c (M°,M 1 ) 

On the other hand, let p = 
Then, we have that 



dip 
in eq. ([T 



1 

4" 



(72) 

where \<f>) = -J=(|01) - |10)). 



distg(M°, M 1 ) > J dip 



h 
4 



1 

2" 



(73) 



D The role of r(B) 

The goal of this Appendix is to prove the following result. 

Theorem 1. Let B be a set of states in B(Hb), and let Hb € L{T-Lb) be a 

non-degenerate energy operator admitting a decomposition in terms of max- 
imal chains j as 



Hb = Y, E (^ + kA)\j,k)(j,k\ 

j k=0 



(74) 



Then, 



e c (M B ) = e Q (M B ) = -{l-r(B)}, 



(75) 



where 
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r(B)=maxJ2 E I0'.* + 1|P|J',*>I- (76) 

p j k=0 

Moreover, for any general one-qubit POVM M = {M x } x , the state p* G B 
maximizing [76] ) allows to generate an effective two-level POVM M = {M x }, 
with 

(M x ) ab {l + 6 ab (l-r(B)}. (77) 

Proof. Let p G B and consider the problem of maximizing 2|(M )io| over all 
M G M{{p}). From eq. (fl9j) . the result is 



2 max 



J 10 



j fc=0 



(78) 



It is clear that the physical POVMs {M jk } maximizing the last expression 
must be chosen such that 



1 ^nh 

1 e 3 K 



=[ ^ I , (79) 

2 2 / 

with e -w ** (j, k + l\p\j, k) eR + U {0}. 

The maximum value of 2|(Mo)io] attainable with B as a resource thus corre- 
sponds to t(B), as defined in eq. ( 1221) . 

Now, let p* G i3 be the minimizer of eq. (122 p [or a very good approxima- 
tion, since the maximum may not be achievable], and suppose that we wish 
to approximate the POVM M G M. Take the physical POVMs to be 

w * = ( I A ) M *(l X ) • ^ 

for all = 0, d. From eq. ffl9l) we have that the final POVM M is then 
given by 



m 

(i\4) 00 = (M x .)oo J2J2 tT ^' fc Xi> fc|) = ( M *)oo, 

j fc=0 
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m 

(M x )n = (M x ) n ]T £ tr(p\j, k - l)(j, k - 1|) = (M x ) n , 

j fc=0 

(M^xo = (M,) 10 ]T ]T e** (j, fc|p|j, fc - 1) 

j fc=0 

= (M x ) 10 t(B), (81) 

and, analogously, (M x )oi = (M x ) 01 r(B). The second part of the theorem is 
thus proven. 

Let us now demonstrate the first part. Consider M = {(I+(— l) a a x ) /2} e 
M.. Then, by the triangle inequality, for any iV G A4(B, 2), we have that 



dist c (M, N) — - max V \ti{p{M x - N x )}\ + V |tr{ P 7V x }| > 
z p * — ' * — ' 

x=0,l x>2 

> imax V |tr{p(M x -^)}| =dist c (M,JV'), (82) 
2 p ^ — ' 

i=0,l 

with N' e M{B) defined as N' Q = N , N[ = J2 x>i N x . 

It follows that, in computing distc(M, A4(B, 2)) we can restrict to two- 
outcome POVMs iV e M(B, 2). For any such POVM we have that 



dist c (M,M(B)) > dist c (M,N) = \ max V \tr{p(M x - N x )}\ = 
= ||M -A^olloo > 



2 p 

21 = 0,1 



i-sup{(iVo)oi :NeM(B, 2)} 



= ^{l-r(B)}. (83) 

Finally, call / : Ai — > M. the map which transforms each POVM M into 
M, as defined in the theorem, and note that distg(M, N) can be written as 

dist Q (M,A0 = l ma xY i tT{p(M x -N x )®S x }, (84) 

with -I < < I. Then, for any POVM M G M, we have that 
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dist Q (M, M{B, 2)) < I max V tr{p(M a , - f{M x )) ® 5"} 

a: 

= maxV|(M a .)oi|tr(p^ ® S*) 

x 

< 1 ~ 2 r(g) ^v / p4Vv / p4)^ 

<^[l-r(B)], (85) 

where i? x = e* 6 *^ 1 0) (1 1 + h.c, for some 9 X G [0, 2ti) and the last inequality 
follows from the fact that the vectors v x = v/ (M X ) Q0 , w x = \J (M x ) n are 
unitary, since Y, x ( M x)oo = J2 x i M x)n = 1- 

From eqs. ( 18 3 p and ( 18 5p . we therefore have that 

\{l - t(B)} < eo(B) < e (B) < i{l - r(5)}. (86) 
The theorem is proven. 

□ 



E Computation of r(C d ) 



It is easy to see that the solution of eq. ( 126|) corresponds to the maximum 
eigenvalue of a d x d matrix of the form 



f 1/2 \ 

1/2 1/2 

1/2 1/2 

1/2 1/2 

V '•• J 



(87) 



Let = Y^k=o c k\k) be an eigenvector of A d with eigenvalue \i. Then, 
the coefficients {c n } must satisfy: 



1 



(c fc+ i + Cfc_i) = /ic fc , for = 0, d - 1, 



(88) 
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with c_i = c<j = 0. Let us try the ansatz c& = sin[(A; + l)a\. Then, condition 
(ISSjl translates as 

-{sin[(/c + 2) at] + sinffca)]} = cos(a) sin[(fc + l)a] = n sin[(fc + l)at], (89) 

with Cd = sin[(d + l)a] = 0. This last condition implies that a = mir/ (d + 1), 
with m G Z. From eq. (l89l) we thus have that Ad has d eigenvalues given by 



Since Ad cannot have more than d eigenvalues, it follows that its whole 
spectrum is contained in ([9"0l . The maximum eigenvalue is obtained by 
taking m — 1, in which case the corresponding (normalized) eigenvector is: 




(90) 




(91) 



F Computation of r(E) 



The aim of this section is to solve the following problem: 



max 




fc=0 



oo 



oo 



S. t. 




(92) 



k=0 



k=0 



where, for simplicity, we have defined E = E/A. 
The dual of this problem is: 



//(A) = — inf 



(V#> = i 



s. t. 



(93) 



with 
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Hb A 

H X = J2 k\k)(k\ -\J2(\k)(k + l\ + \k + l)(k\). (94) 

k k 

For finite values of A, the infimum in eq. f[9"3"j) is actually a minimum, 
i.e., problem (|93|) can be reformulated as an eigenvalue problem (namely, 
computing the minimum eigenvalue of H\). To see this, note that we can 
approximate the above minimization over infinite-dimensional vectors \ip) = 
SfcLo c k\k) by finite dimensional optimizations over vectors of the form \%p') G 
C d . Indeed, let |^) = Y^=o c k\k) be a normalized vector with (ip\H\\i(:) < 
(it is easy to see that //(A) > for all A > 0). Since = 2, we have that 

< -(i/>\H x \iJ>) < -(ip\H B \^) + 2A. (95) 
It follows that J2 n \ c n\ 2n < 2 A, and so, for any DeN, 

^|c„| 2 <^|c n | 2 n<-^ 1 . (96) 

n>D n>D 

Under the assumption that c n > for all n, it is easy to see that 



D+l 

Wm _ g < (Vto|A|Vto) < (97) 

D+l 

where I^d) oc X^n<D c "l n ) e C D+1 is a (normalized) finite dimensional ap- 
proximation of Since we can make D arbitrarily large, it is obvious that 
problem (1931) can be approximated arbitrarily well by its finite-dimensional 
analog 



min (^\Hf\ip) 
s. t. |V) G C d , (V#> = 1, (98) 

with H d x = Eto k\k)(k\ - A Etl(\k)(k + 1| + |*><* + 1|). 

Call the minimizer of this latter problem, and consider the sequence 
of minimizers (|V>d))d- Then, there exists a subsequence (di,c?2, •••) such that 
lim^->oo(0|^.) = Cq converges. Likewise we can find a subsequence (eQj of 
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(cQj such that lim^_ >0O (l|^i,) = = c\ also converges. Iterating this 

1 i 

process, we obtain a sequence of values (cq, c*, ...). Now, as (ipd\4>d) = 1 
and the weight of the terms ^ n>Z) n|(n|^)| 2 becomes negligible for large 
.D irrespective of cf, we have that the vector \ip*) = X]fcLo c fcl^) satisfies 
|||^)|| 2 = land(^|ifA|^) = /i(A). 

Let us thus solve f )98|) . Let \ip) = J2k=o Ck \^) ^ e an eigenvector of i?^ 
with eigenvalue — /x, with /x > 0. Then, 



k + fi 

Cfc+i = - c fc - c fc _i, (99) 

with c_i = q = 0. 

Now, for A > 0, it is clear that the eigenvector with minimum eigenvalue 
of Hf can be chosen such that Ck > 0. It is also straightforward that it 
cannot contain intermediate zeros in its component vector, that is, it cannot 
be of the form 

K d-l 
k=0 k=K+L 

for any L > 1. Indeed, if such were the case, the vector 

K d-l 

W) =Y, c k\k) + Ck\k-L + l) (101) 

fe=0 k=K+L 

would satisfy (^'\Hf\ip') < (ip\H%\ip). Finally, by the recurrence condition 
f )99|) . it can also be seen not to be of the form (ci, C2, ck, 0, 0, 0) for 
A > 0. It follows that Ck > for = 0, d — 1. For = 0, d — 1 we can 
thus write the c's as = niLo 3 -*' f° r some Xk > 0. Substituting in eq. ( l9"9"j) 
we have that 

^ = Ma 1 r ' (102) 

for k > 1. We also derive the conditions Xi = ^, = 0. 
Iterating (11021) . we have that 

f = *i = i + „ 1 | — (103) 

A ~~A 2 + m 1 
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Since = 0, this implies that the minimum eigenvalue /x of Hf mus satisfy 
the characteristic equation 

^ — ^ = 0. (104) 



A 1±M 



Moreover, as long as any // satisfies eq. (I104p with each term in the 
fraction being different from zero, we can construct a set of weights 
such that (l,xi,x\X2,xiX2X3, ...) is an eigenvector of Hf with eigenvalue /i. 
Suppose now that \i fulfills eq. (I104p . but some intermediate subfraction 



becomes zero, i.e., 



k + \i 



fc+l+M 



0, (105) 



d-l+n 



with d — 1 > k > 1, and such that no intermediate subfraction of (11051) is 
zero. We could then create a vector (1, xi, Xix 2 , n^fc 1 x i) e C d ~ fe which 
happens to be an eigenvector of H^~ k + kl with eigenvalue —ji. But the 
minimum eigenvalue of H^~ k + kl is strictly greater than that of Hf. 
From the above discussion, it follows that the minimum eigenvalue of Hf is 
minus the maximum fi that satisfies eq. (11041) . 

What happens when we take the limit d — > oo in expression (11041) ? In 
[23] . page 349, we find the following beautiful identity: 

J n-i{z) _ 2n _ 2(^+1) n nfi x 

'* V / 1 (n + l)(n + 2) 

^ — 



1 (n + l)(n + 2) 
1- ' • 

Multiplying each numerator and denominator by 2 ("+ fc ) accordingly, the 
latter expression can be shown equal to 

2n 1 , s 

T " ^ i — (107) 

z 2(n + 2) 1 

z 2(n + 3) j~ 
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Identifying A = z/2 in (11041) . we have that the solution of (I9"3"j) is the 
maximum /i > such that 



(108) 

By the interlacing properties of the zeros of the Bessel function, this last 
condition is equivalent to J M _i(2A) = 0, and so the solution of problem ( 193|) 
is given by the solution of the transcendent equation 

= 2A. (109) 

From now on, we will denote its solution /i(A). 

How does this relate to our initial problem ( 192]) ? Let be such that 
(ip\Hg\ip) = E. By definition, we have that 

H x + fi(\)>0, (110) 

for all A > 0. 

Bracketing this expression by we conclude that 

\CkCk+l\ < • (HI) 

k 

Minimizing with respect to A, we thus arrive at 

J2\ckC k +i\<^{E), (112) 

k 

with (f(E) defined as in (13"U1) . To show that this last inequality is tight, denote 
by the eigenvector with minimum eigenvalue of H\ and notice that 



- fi(X + 6X) + /i(A) = — — - {(^a+<5a|#a+«5a|^a> - (^x+sx\H x \^x)} 

\yx+6x\yx) 

Dividing by 5X and taking the limit 5X — > 0, we have that 

MM^Px) = |^(A). (114) 
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It follows that 



£?(A) = ^aI^a) = A^-MA). (H5) 

Now, it is easy to see that E(0) = and that E(X) is continuous in A. 
On the other hand, we have that, for jtt> 1, 

j lt , 1 =H + c^ 3 + 0(n- 1 ^ 3 ), (116) 

with c ps 1.85575 [2i] . The identity = 2A thus implies that /i ps 2A — 

2 i/3 cA i/3 + Substituting in expression f)115p . we have that (i()\\Hb\'4 } \) 

diverges for A — > oo. All this implies that, for any E > 0, we can find A > 
such that {ip\\H B \ip\) = E. Such a state clearly saturates inequality fillip . 
The relation r({p : tr(pife) < -£/}) = <p(E) is thus proven. 
It is worth noticing that any value of A = A locally minimizing the right hand 
side of f illip must satisfy 

a * + "«=0. (117) 



d\ 2A 

Or, equivalently, E = Xdfi(X)/dX — /x(A). The state \ip\) thus saturates 
inequality f illip for A = A, and, consequently, for A ^ A the r.h.s. of fillip 
cannot become smaller. Any extreme point of the latter function is hence a 
global minimum, and so computing <f>(E) numerically becomes an easy task. 

Finally, let us speak about the asymptotic behavior of <p{E). Replacing 
/i(A) by 2A — 2 1 / 3 cA 1 / 3 in the right hand side of eq. f illip , and minimizing 
with respect to A, we find that A* ~ And, consequently, 

<^)«1-^- (H8) 



Substituting the value of c, we arrive at (J3T 
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